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ABSTRACT 



The properties of the stellar populations in ULIRGs II: 
star formation histories and evolution 
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o 

' This is the second of two papers presenting a detailed long-slit spectroscopic study 

of the stellar populations in a sample of 36 ULIRGs. In the previous paper we presented 
i-^ ■ the sample, the data and the spectral synthesis modelling. In this paper we carry out 

a more detailed analysis of the modelling results, with the aim of investigating the 
q . general properties of the stellar populations (i.e. age, reddening and percentage contri- 

' bution) and the evolution of the host galaxies, comparing the results with other studies 

c/3 \ of ULIRGs and star forming galaxies in the high-z Universe. The characteristic age of 

C$ . the young stellar populations (YSPs) is < 100 Myr in the nuclei of the overwhelming 

majority of galaxies, consistent with the characteristic timescale of the major burst 
of star formation associated with the final stages of major galaxy mergers. However, 
the modelling results clearly reveal that the star formation histories of ULIRGs are 
' complex, with at least two epochs of star formation activity. Overall, these results are 

, consistent with models that predict an epoch of enhanced star formation coinciding 

■ with the first pass of the merging nuclei, along with a further, more intense, episode 

of star formation occurring as the nuclei finally merge together. It is also found that, 
although YSPs make a major contribution to the optical emission in most of the ex- 
tended and nuclear apertures examined, they tend to be younger and more reddened 

Q\ \ in the nuclear regions of the galaxies. This is in good agreement with the merger sim- 

■ ulations, which predict that the bulk of the star formation activity in the final stages 

of mergers will occur in the nuclear regions of the merging galaxies. In addition, our 
results show that ULIRGs have total stellar masses that are similar to, or smaller than, 
the break of the galaxy mass function (i.e. ULIRGs are sub-m» or ~ m* systems), 
and that the young stellar populations detected at optical wavelengths dominate the 
stellar mass contents of the galaxies. Finally, we find no significant differences between 
the ages of the YSP in ULIRGs with and without optically detected Seyfert nuclei, 
nor between those with warm and cool mid- to far-IR colours. While this results do 
not entirely rule out the idea that cool ULIRGs with HII/LINER spectra evolve into 
warm ULIRGs with Seyfert-like spectra, it is clear that the AGN activity in local 
Seyfert-like ULIRGs has not been triggered a substantial period (>100 Myr) after the 
major merger-induced starbursts in the nuclear regions. 
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1 INTRODUCTION 

First discovered in large numbers by the IRAS satellite, ul- 
traluminous infrared galaxies (ULIRGs, Lir, > 10 12 Lq) rep- 
resent some of the most rapidly evolving objects in the local 
Universe. These objects are important in several respects, 
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such as testing the merger models, investigating the nature 
of the AGN-starburst connection in merging systems, and 
understanding the physical processes taking place in the re- 
cently discovered high-redshift star forming galaxies, which 
have similar properties to the nearby ULIRGs. 

A remarkable property of ULIRGs is that more than 
90% of them are associated with galax y mergers and in- 
teractions fsee lSanders fe Mirabell . Il996l . for a review). The 
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simulations have shown that, during galaxy collisions in- 
volving gas-rich disk galaxies, the gas loses angular momen- 
tum due to dynamical friction, tidal torques, and inflows 
towa rds the centres of the galaxies (e.g. Mi hos fc Hernquistl 
1996). Such high nuclear concentrations of gas are ca- 
pable of triggering both starburst and AGN activity. In 
general, these simulations predict two epochs of enhanced 
star formation in major galaxy mergers: the first occurring 
around the first encounter, and t he second episode when the 
nuclei are close to c oalescence jMihos fc Hernquistl 199€ 
Barnes fc Hernquistl. 119961; jSpringel et all 120051 ; I Cox et al. 



effective radii and velocity dispersions smaller than those 
of their comparison sample of nearby radio galaxies. Fur- 
ther more, the dynamical masses of UL I RGs are 10 10 - 10 11 
Mr^ dTacconi et all |2002| ; IColina et all |2005| ; iDasvra et ail 



20061 ; Idi Matteo et all 120071 ) . However, both the time lag 



and the relative intensity of the peaks of starburst activ- 
ity during the merger event depend on several factors, such 
as the presence of bulges, feedback effects, gas content and 
orbital geometry. Given that the models make specific pre- 
dictions about the histories of the star formation triggered 
in the course of major gas-rich mergers, studies of the stel- 
lar populations in ULIRGs provide useful information about 
the mergers and, potentially, allow us to test the models. 

ULIRGs are classified on the basis of their infrared 
colours as warm (/25//60 > 0.2 J] and cool (/25//6O < 0.2) 
ULIRGs. Warm ULIRGs represent ~20 - 25% of the total 
population of ULIRGs discovered by IRAS. Most of such 
objects have an AGN optical spectral classfication, tend to 
be more compact than cool ULIRGs , and are frequentl y 
found in an advanced merger state l|Surace et all 1 1998T) . 
These properties suggest that ULIRGs play an important 
role in the formation and evo lution of QSOs (Sande rs et 
al., 1 988b) and radio galaxie s l|Tadhunter et all l2005h . In- 
deed, ISanders et al. (1988b)l proposed a scenario in which 
cool ULIRGs evolve into warm ULIRGs on their way 
to becoming typical optically-selected QSOs. Further evi- 
dence in f avour of such an ev o lutionary scenario was re - 
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(Il998ll. ISurace fc Sanders! (|l999h . 
ISurace et all i|200d) and lSurace fc Sanders! (|2000h on the ba- 
sis of their imaging s t udies o f c ool an d warm ULIRGs, and 
ICanalizo fc Stocktonl (|2000al lbl, l200ll ). based on a spectro- 
scopic analysis of a sample of transition QSOs — objects 
that, based on their infrared colors, may represent a tran- 
sitionary stage between ULIRGs and QSOs. More recently, 
iTadhunter et alj (|2007T l carried out a mid- to far-IR study 
of powerful radio galaxies, searching for signs of hidden star 
formation activity. They found that the fraction of powerful 
radio galaxies with energetically significant star formation 
activity (20% - 30%) is similar to that deduced from opti- 
cal studies. The derive d ages and masses of the young stel- 
lar po pulat ions in the Wills ct al. (|2004T ). ITadhunter et al.l 
(2005) and ITadhunter et al.l 1 2007 ) samples of radio galax- 
ies, and the fact that some of these radio galaxies are actu- 
ally classified as ULIRGs, is consistent with the idea that 
some radio galaxies ar e the evolved product of ULIRGs 
ITadhunter et all 120051 ). 

In terms of the end products of the merge r events asso- 
ciated with ULIRGs, the kinem atic studies of iGenzel et alj 
l|200ll ) and ITacconi et all (|2002h have shown that ULIRGs 
match the locations of intermediate mass, disky ellipticals 
and SO galaxies in the fundamental plane. ULIRGs have 
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which is similar to intermediate-mass elliptical 
galaxies. Based on the dynamical properties of ULIRGs, 
ITacconi et ail l|2002l ) concluded that not all ULIRGs can 
evolve into typical QS Os or powerfu l radio galaxies. How- 
ever, the recent work of lDasvra et al.l (|2007T ). based on near- 
IR spectroscopic observations of a sample of 12 (mainly 
Palomar Green) QSOs, has shown that it is possible for 
ULIRGs to have similar dynamical properties to QSOs. 
Therefore, it is clear that, despite the recognition of an evo- 
lutionary link between ULIRGs, QSOs/radio galaxies and 
elliptical galaxies, considerable uncertainties remain about 
the true nature of the link. 

Studies of ULIRGs are also important in the context 
of the various populations of high-z galaxies that have 
been re cently detect e d usi ng near-IR (distant red galaxies, 



DRGs, iFranx et all l2003h. mid-IR (24 / zm Spiteer-selected 
galax ies . ICaputi et al. , 2006; lYan etafl 120071 ISaiina et all 
120071 . and others), a nd sub-millimeter wavelengths (sub- 
mm galaxies, SMGs, ISmail et all |2002| ; iBlain et all |2002| ; 
IChapman etail 120031 ; iPope et all |2005| . and others). These 
objects have properties similar to the nearby ULIRGs (stel- 
lar mass and luminosity, for example), and make up a sub- 
stantial fraction of the star formation (SF) at z ~ 1-2. As 
local analogues of such objects, ULIRGs provide an oppor- 
tunity to study the physical processes associated with their 
prodigious SF activity in depth. 

It is also important to emphasize that most of the stud- 
ies of the stellar populations in high redshift galaxies are 
based on the modelling of a few photometric points and, 
therefore, it is hard to constrain the properties of these pop- 
ulations. Sensitive studies of star forming galaxies at lower 
redshifts are useful, not only to better understand the links 
between local starburst objects and their high-z counter- 
parts, but also to gain valuable information on the uncer- 
tainties associated with the derived properties of the stellar 
populations in high-z star forming galaxies. 

At this stage it it is clear that there are many unresolved 
issues concerning the nature of ULIRGs. In particular: 

• Do the merger simulations adequately describe the ob- 
served properties of ULIRGs? 

• Do the properties of the stellar populations in ULIRGs 
correlate with other properties of ULIRGs, such as infrared 
luminosity, interaction class or spectral classification? 

• What is the nature of the links between ULIRGs and 
other types of object in the local universe, such as QSOs, 
radio galaxies and elliptical galaxies? 

Given their importance for studying the evolution of 
galaxies via major gar-rich mergers, it is perhaps surpris- 
ing that there have been relatively few studies of stellar 
populations in ULIRGs. Many previous studies have con- 
centrated on studying the SEDs and mid-IR emission line 
spectra of ULIRGs, with an emphasis on determining the 
dominant heating source for the MFIR-emitting dust (e.g. 
iFarrah et all [2003 ). However, such studies are not capable 
of determining the detailed properties of the stellar pop- 
ulations and how they vary with spatial location in these 
complex systems. Therefore, order to address the key issues 
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surrounding the evolution of ULIRGs we have undertaken 
a programme of deep optical long-slit spectroscopic obser- 
vations of a substantial sample of nearby ULIRGs, aimed 
at investigating their star formation histories. In Rodriguez 
Zaurm et al. (2009a), hereafter Paper I, we presented the 
sample, data reduction and spectral syntehis modelling tech- 
niques. In this paper we discuss the results presented in 
Paper I in the context of evolutionary models for merging 
systems and the results obtained in studies of stellar popu- 
lations in high-z star forming galaxies. 

Throughout this paper, we assume a cosmology with Ho 
= 71 km s" 1 Mpc~\ D.0 = 0.27, Q A = 0.73. 



2 MODELLING THE STELLAR 

POPULATIONS: GENERAL RESULTS 

In Paper I we presented the results of spectral synthesis 
modelling of the complete (CS) and the extended (ES) sam- 
ples, including 26 and 36 ULIRGs respectively. The de- 
tailed properties of these samples are described in that pa- 
per. To summarize, in the first place, we observed a com- 
plete RA-limited and declination-limited sub-sample of the 
iKim fc Sandersl (|l998h 1 Jy sample of ULIRGs, with RAs in 
the range 12 < RA < 1:30 h, declinations 5 > -23 degrees, 
and redshifts z < 0.13. The redshift limit was chosen to en- 
sure that the objects are sufficiently bright for spectroscopic 
study, and also to keep the size of the sample tractable for 
deep observations on 4m class telescopes. We will refer to 
this sample of 26 objects as the complete sample (CS, see 
Table 1 in Paper I for details of some of the the properties 
of the ULIRGs in the CS). In addition, we also observed 10 
objects outside the redshift range and/or the RA and Dec 
range of the CS. With the aim of including a larger num- 
ber of warm objects (only 8 warm ULIRGs are included in 
the CS), 7 of these 10 ULIRGs have warm mid- to far-IR 
colours. The sample including the CS and these 10 addi- 
tional objects will be referred to as the extended sample 
(ES). Thus, the ES comprises 36 objects: 21 cool ULIRGs, 
and 15 warm ULIRGs. Since the ES accounts better for the 
diversity within the ULIRG phenomenon, and our statis- 
tical analysis shows no significant differences between the 
properties of the stellar populations in the ES and CS, we 
concentrate on the ES sample for the study presented in this 
paper. 

In order to perform a detailed study of the stellar pop- 
ulations in ULIRGs, a total of 133 apertures was extracted 
for the objects in the ES. The extraction apertures were se- 
lected from spatial cuts of the 2-D frames in the line-free 
continuum wavelength range 4400 - 4600 A, based on the 
visible extended structures and the requirement that the 
apertures are large enough to have a sufficiently high S/N 
ratio for further analysis. To compare the stellar populations 
between the objects in our sample, apertures with a metrical 
scale of 5kpc centred on the main nuclei were extracted for 
all the objects in the extended sample (ES), including sepa- 
rate extractions for multiple nuclei in individual sources. A 
second set of apertures was then selected to sample the spa- 
tial features of those objects in the ES showing tails, bridges 
and other diffuse structures. 

To perform the fits, we have used the CONFIT code 
(see lRobinson et al.ll200ol and RZ08 for details). The CON- 



FIT approach consists of a direct fit of the overall contin- 
uum shape of the extracted s pectra using on a minimum \ 2 
techn ique l)Tadhunter et all 120051 : iRodrfguez Zaurm et all 
2007). CONFIT is based on a "simplest model" approach, 
i.e. we fit the minimum number of stellar components re- 
quired to adequately model the data. Therefore, CONFIT 
allows for a maximum of two stellar components plus a 
power-law in some cases. Throughout this paper, we de- 
fine young stellar populations (YSPs) as stellar compo- 
nents with ages tysp< 2 Gyr, and old stellar populations 
(OSP) as components with ages tosp>2 Gyr. With the aim 
of better describing each combination, it is convenient to 
further sub-divide the YSPs into two groups: very young 
stellar populations (VYSP): stellar components with ages 
twsp < 0.1 Gyr; intermediate-age young stellar populations 
(IYSP): stellar components with ages in the range of 0.1 < 
tiYSP < 2 Gyr. 

Given that the properties of the stellar populations were 
not know a priori, three modelling combinations were inves- 
tigated. The results obtained for each of the combinations 
are described in detailed in Paper I, and are summarized as 
follows. 

• Combination I. As a first approach we used a two 
component model comprising a YSP (tysp < 2 Gyr) with 
varying reddening (0.0 < E(B — V) < 2.0, increasing in steps 
of 0.1), along with an unreddened OSP of age 12.5 Gyr. 
Accounting for an old component covers the case in which 
one or more of the merging galaxies is early-type or bulge- 
dominated galaxy. Using this combination we obtained ad- 
equate fits (Xrcd < 1-0) to the overall shape of the SED 
for all but one (IRAS 23060+0505, for which a power-law 
is required) of the extraction apertures modelled. However, 
Comb I fails to adequately fit the absoption detailed absorp- 
tion - particularly the Call K feature - for ~ 42% (54) of 
the apertures modeller^. 

• Combination II. This combination consists of three 
components: an OSP of age 12.5 Gyr and zero reddening, 
along with a YSP with variable reddening (0.0 < E(B— V) < 
2.0, increasing in steps of 0.1), and a power-law with a spec- 
tral index in the range —15 < a < 15. The power-law is 
included to repres ent either scatter e d or d irect AGN contin- 
uum component l|Tadhunter et all |2002| ). or a highly red- 
dened VYSP. Adequate fits (y 2 < 1.0) to the overall SED 
shapes were found for all of the apertures of all objects us- 
ing this combination. On the other hand, this combination 
fails to fit the detailed absorption features for 10% (13 aper- 
tures) of the apertures modelled. We find that the minimum 
precentage contribution of the OSP component is less that 
10% (due to the uncertainties inherent in the modelling tech- 
niques percentage contributions < 10% are considered neg- 
ligible) for most of the apertures. This result suggest that 
an OSP component is not required to fit the optical spectra 
of the majority of the objects in the ES. 

• Combination III. A drawback of Comb II is that, for 
those cases in which the power-law is likely to represent a 
VYSP component, it gives no information about the detailed 
properties (age and reddening) of such a component. In ad- 
dition, for ULIRGs as a class, there is no reason, a priori, to 

2 As discussed in paper I, some of this discrepancy may be due 
to the effect of ISM absorption on the Call K resonance line. 
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expect a major contribution to the optical emission from a 
12.5 Gyr OSP. Such lack of an OSP is also suggested by the 
results of Combination II. Therefore, in order to explore the 
possibility of YSPs dominating the optical emission from 
the objects we decided to use a combination consisting of 
two YSP components: a IYSP with ages in the range 0.3 
- 2.0 Gyr, with E(B - V) values of 0.0, 0.2 or 0.4, plus a 
VYSP with age in the range 1 - 100 Myr and reddenings 
(0.0 < E(B-V) < 2.0, increasing in steps of 0.1). Adequate 
fits to both the continuum and the detailed absorption fea- 
tures are obtained for all but 9 (~ 7%) of the extraction 
apertures modelled. In the minority of the cases for which 
no adequate fits could be obtained using Comb III, either 
an OSP or a power-law component was required in addition 
to a YSP. 

Signifcantly, for the majority of apertures we find that 
both the SEDs and the detailed absorprion features can be 
adequately modelled with all three combinations, although 
in most cases Comb III provides the best overall fit. This 
demonstrates that it can be difficult to derive a unique model 
solution even with high-quality, wide spectral coverage opti- 
cal spectra. Nonetheless, the fits provide the following gen- 
eral results. 

- Dominant young stellar populations. YSPs dom- 
inate the optical light, while OSPs are relatively unim- 
portant, in the overwhelming majority of apertures mod- 
elled. In terms of the nuclear apertures, the only excep- 
tions are IRAS 13451+1232 5kpc, IRAS 1648NE 5kpc, 
IRAS 21208-0519 5kpcII, IRAS 23233+2817 5kpc and 
IRAS 23327+2913 5kpcII. Although there is some ambigu- 
ity in the cases of IRAS 13451+1232 5kpc, IRAS 1648NE 
5kpc, and IRAS 23233+2817 5kpc which can be modelled 
either in terms of a dominant OSP plus a YSP (Comb I) or 
a dominant "old" IYSP (>1 Gyr) plus a VYSP (Comb III), 
we can only adequately model IRAS 21208-0519 5kpcII and 
IRAS 23327+2913 5kpcII with Comb I models that include 
an OSP. Note that, in the case of IRAS 23327+2913 5kpcII, 
it is possible to model the optical spectrum using an OSP, 
without any requirement for a YSP. Apart from these ex- 
ceptions, for most objects the results are consistent with the 
idea that the precursor galaxies are gas-rich spiral galaxies. 

- Complex star formation histories. Consistent 
with our previous resul t s for the nearby ULIRG Arp220 
(Rodrigue z Zaurin et all I2008T ), we find that most aper- 
tures of most objects are best fitted with models that in- 
clude a combination of VYSP and IYSP (i.e. Comb III). 
This suggests that the recent star formation histories of 
most ULIRGs are complex, with at least two star formation 
epochs, even if the most recent (represented by the VYSPs) 
often dominates the optical light. 

In the following sections we consider the properties of 
the YSP in ULIRGs in more quantitative detail, and inves- 
tigate whether they are correlated with other properties of 
these merging systems. 

2.1 Characteristic properties of the YSPs in the 
nuclear and extended regions. 

As described in Paper I, Comb I (YSP + OSP) provides use- 
ful estimates of the "luminosity weighted" ages, reddenings 
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Table 1. Mean and median values for the parameters associated 
with the YSPs for both the nuclear and extended regions of the 
objects in the ES, obtained using Comb I. 



D Confidence 
Level 

LW-age 0.325 99% 

LW-E(B - V) 0.50 > 99.9% 
LW-% 0.383 > 99.9% 

Table 2. Results of the K-S test for the significance of differ- 
ences between the distributions in the nuclear and the extended 
regions, of the LW-age, LW-E(B - V) and LW-% for the ES and 
using Comb I. Col (2): the maximum deviation between the two 
cumulative distribution functions. Col (3): the probability that 
the difference between the distribution of the nuclear and the 
extended regions is significant. 

and percentage contributions of the dominant YSPs within 
each galaxy. Providing a single estimate (with an associated 
uncertainty) for the parameters associated with the YSPs 
in each aperture, this combination is particularly useful for 
performing a statistical analysis of the properties of the stel- 
lar populations in the CS and the ES samples. 

Figure [l] shows the distributions of the average 
luminosity- weighted ages (LW-age), reddenings (LW-E(B — 
V)) and percentage contributions to the flux in the normal- 
ising bin (LW-%) of the YSPs in the nuclear and extended 
apertures. This figure has been made using 34 of the 36 ob- 
jects in the ES. For two of the three objects classified as Syl 
(IRAS 15462-0450 and IRAS 21219-1757) no extended aper- 
tures were extracted and, therefore, they are not included 
in the analysis presented here. In the case that one of the 
extended apertures samples a similar region to the 5 kpc 
aperture, only the latter is used for the plots (for example, 
ApE and 5 kpc in the case of Mrk 273). In the case of Arp 
220, only two of the 24 apertures were included, in order to 
avoid giving too much weight to the extended apertures in a 
single object. These are the 5 kpc aperture and APtotal for 
PA 75*. Overall, a total of 100 of the 128 apertures modelled 
for the ES were eventually used for the figures. 

The upper panel of Figure \T\ shows the presence of a 
clear upper limit of 100 Myr for the LW-ages of the YSPs 
in the nuclear apertures; only ~ 9% of the nuclear apertures 
(4 of the 44 nuclear apertures used) have average LW-ages 
older than this value. In contrast, ~41% of the extended 
apertures (23 of the 56 extended apertures used) have av- 
erage LW-age values older than 100 Myr, and as old as 1.2 
Gyr. In terms of reddening, Figure [1] (middle) suggests the 
presence of higher reddening (LW--E(B — V)) values in the 
nuclear regions of the galaxies compared with the extended 
regions. We also find that the YSPs represent the dominant 
contribution (LW-% > 50%) to the optical emission for all 



Optical spectroscopy of ULIRGs: star formation histories and evolution. 



5 



44- 
40- 
36- 
32- 
28- 

| 20- 
Z 16- 
12- 



NUC AP 



-Q- 



400 600 800 1000 
LW Age (Myr) 



44 
40 
36- 

32 

28 

55 24 
-o 

S 20 

=3 

Z 16- 

12 
8 
4 




EXT AP 



200 400 600 800 1000 
LW Age (Myr) 



20-p 
18- 
16- 
14- 
12- 
I 10- 

s 

I 8 : 

6- 

4-- 



NUC AP 



0.4 0.6 0.8 1.0 1.2 

LW-E(B - V) 



32 

28 

24 

20 

I 16 
6 

Z 12-1 

8 
4 



0.4 0.6 0.8 
LW-E(B - V) 



EXT AP 




10 20 30 40 50 60 70 
LW% 



12- 
10- 



EXTAP 



10 20 30 40 50 60 70 
LW% 



Figure 1. Histograms showing the LW-age (top), LW-E(B - V) (middle), and LW-% (bottom) contribution distributions for the nuclear 
(left) 5 kpc and the extended (right) apertures. The results of Comb I show that YSPs are present at all locations in the galaxies for the 
majority of the objects in the ES, with a clear cut-off in YSP ages at 100 Myr in the case of the nuclear regions. The figure also suggests 
that the YSPs in the nuclear regions of the galaxies, coinciding with higher reddening, represent a more important contribution to the 
flux at optical wavelengths than those in the extended regions. 



but 11 (11%) of the apertures used (4 nuclear apertures and 
7 extended apertures, Figure fTJ bottom panel). However, 
the YSP make a more important contribution to the optical 
light in the nuclear regions. 

Table [1] shows the mean and median values of the dif- 
ferent parameters associated with the YSPs. To investigate 
whether the distributions of the YSPs properties in the nu- 
clear and extended regions are indeed different, we used the 
Kolmogorov-Smirnov (K-S) two sample test, and the results 



are presented in Table [2] We find that the differences be- 
tween the LW-age, LW-E(B - V) and LW-% distributions 
for the YSPs in the nuclear and extended regions of the 
galaxies are significant with a confidence level of > 99% (i.e. 
the probability that the distributions are drawn at random 
from the same parent distribution is <1%). 

To further investigate the properties of the stellar pop- 
ulations, we performed a similar study to that presented 
above using the results of Comb III (VYSP + IYSP). Fig- 
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Figure 2. Histograms showing the distributions of the average VYSP- Age (top), E(B - V) values (middle) and the minimum percentage 
contribution to the flux in the normalising bin (Pvysp > bottom) for the nuclear 5 kpc (left) and the extended (right) apertures, using 
the ES and Comb III. The locations of Arp 220 and PKS1345+12 are also shown in the figure. 
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Table 3. Mean and median values for the parameters associated 
with the VYSPs, for both the nuclear and extended regions of the 
objects in the ES, obtained using Comb III. 



ure [2] shows the distributions of the average values for the 
VYSPs age, reddening and the minimum percentage con- 
tribution to the flux in the normalising bin (Pvysp )i f° r 
both the nuclear and the extended regions. From a visual 
inspection of the histograms, the modelling results suggest 
the presence of "younger" and more reddened VYSPs in the 
nuclear regions of the objects. Moreover, the lower panel in 
Figure [2] is designed to show the importance of the VYSPs 
within the different regions of the galaxies sampled by the 
apertures. The figure suggests that VYSPs are more impor- 
tant in the nuclear regions of the galaxies. Table [3] presents 
the mean and median values for the VYSP-Age, VYSP-E(B 
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D Confidence 
Level 

VYSP-Age (Myr) 0.27 95% 
VYSP-E(B - V) 0.27 95% 
P^? SP 0.35 99.5% 

Table 4. Results of the K-S test for the significance of differences 
between the distributions in the nuclear and the extended regions, 
of the VYSP-Age, VYSP-E(B - V) and P m ^? S p for the ES, using 
Comb III. Col (1): the maximum deviation between the cumula- 
tive distribution functions Col (2): the probability that the dif- 
ference between the distribution of the nuclear and the extended 
regions is significant. 

- V) and Pvysp, wh ile Table 0] shows the results of the K-S 
test obtained using Comb III. We find that the difference 
between the VYSPs age and E(B - V) distributions in the 
nuclear and the extended regions is significant with a con- 
fidence level of 95%. In the case of Pvysp> the confidence 
level is 99.5%. Overall, the results are consistent with those 
of Comb I, but the level of significance is lower. 

The lower panel of Figures [2] also shows the loca- 
tions of IRAS 13451+1232 (PKS1345+12) and Arp 220, 
the two objects for which deta iled s tudies were presented in 
iRodriguez Zaurin et al.l (|2007t ) and iRodriguez Zaurin et al.l 
12008). Arp 220 is frequently used as the archetype of 
ULIRGs as a class. However, the modelling results indicate 
that the YSPs in this object are not typical of the objects 
ULIRGs in general. While significant VYSPs are detected at 
all locations of the galaxies sampled by the apertures, in the 
case of Arp 220 the contribution of such populations is rela- 
tively small, apart from in the nuclear region of the galaxy. 
Furthermore, in this object, a uniform IYSP of age 0.5 - 0.9 
Gyr is found at all locations sampled by the apertures. 

In the case of PKS 1345+12, some peculiarities might 
be expected, since this is the only object in our sample clas- 
sified as radio galaxy. For this object, either OSPs or dom- 
inant "old" IYSPs (1-2 Gyr) are required to model the 
data (see Paper I for details); the contribution from VYSPs 
is relatively small. Independent of the combination of stel- 
lar populations assumed, this galaxy is the most massive 
amongst the objects in our sample (see Table 4 in Paper I). 

Overall, our modelling results indicate that the VYSPs 
in the nuclear regions of the galaxies make a more significant 
contribution to the optical light than those of the extended 
regions. In terms of age and reddening, the VYSPs located in 
the nuclear regions tend to be younger and more reddened, 
although further analysis is required to confirm these results. 
This results will be discussed in the context of the merger 
simulations in Section 6. 



2.2 Correlations between the properties of the 
YSP and other ULIRG properties 

Given that certain of the general properties of ULIRGs (e.g. 
IR luminosities, mid- to far-IR colours, optical morphology 
and spectral class) may change as the systems evolve along 
the merger sequence, it is interesting to determine whether 
such properties correlate with those of the YSP determined 
from spectral synthesis modelling. 

With the aim of investigating whether the more lu- 
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Table 5. Mean and median values of the VYSP-Age, VYSP-E(B 
-V) and P^ SP for the HII, LINER and Sy2-ULIRGs, for the 
nuclear 5 kpc apertures. 



minous objects have younger, redder or more important 
VYSPs, Figure [3] shows the average VYSP age, reddening 
(E(B - V)) and minimum percentage contribution (Pvysp) 
obtained for the nuclear 5 kpc apertures, plotted against the 
log of the IR-luminosity (from Kim & Sanders 1998) . Note 
that, since it was not possible to model the Syl-ULIRGs 
(IRAS 12540+5708, IRAS 15462-0450 and IRAS 21219- 
1757, see Paper I for details), no Syl galaxies are shown 
in the figure. We find that the properties of the VYSPs are 
independent of the IR-luminosities of the sources. 

In addition to looking for trends with luminosity we can 
also examine whether the stellar properties depend on the 
optical spectral classification. As can be seen from Figure 
[3] where the points are labelled according to spectral clas- 
sification, we see no clear trend with optical type. Table [5] 
shows the mean and median values of the VYSP age, red- 
dening and Pvysp f° r the different types of object i.e. HII-, 
LINER- and Sy2-ULIRGs. If, for example, there exists an 
evolutionary link between the HII- or LINER-ULIRGs and 
the Sy2-ULIRGs (Sanders et al, 1988b), the latter represent- 
ing a more advanced stage in an evolutionary sequence, one 
might expect the Sy2-ULIRGs to have older ages. No trends 
are observed in Figure [3] for the different types of ULIRGs. 
A first interpretation is that there is no evolutionary link 
between HII-, LINER- and Sy2-ULIRGs. However, it is also 
possible that the transition between one type and the other 
occurs over a shorter timescale than the typical timescale 
of the enhancement of the star formation activity in the fi- 
nal s tages of the merger event as the nuclei coalesce (~ 100 
Myr, I Barnes fc Hernquistt ll996l L In that case, no trend is 
expected to be observed. 

We have also investigated the presence of correlations 
between the properties of the VYSPs and other properties of 
ULIRGs such as the Veilleux et al. (2002) interaction class, 
or the MFIR colour ratio /2s //6C0- The conclusions reached 
are identical, i.e. no evidence for correlations or trends is 
found. 

As well as the explanation already discussed above, pos- 
sible reasons for the lack of correlations between the YSP 
properties and other properties of ULIRGs include the fol- 
lowing: 



3 The quantities /2s and /go represent the IRAS flux densities in 
Jy at 25/xm and 60/^m. 
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Figure 3. Average VYSP age, E(B - V), and P^? sp plotted 
against Log(Lip),). The nuclear 5 kpc apertures have been used 
for these plots. The points are labelled based on the Veilleux et 
al. (1999) spectral classifications. No significant correlations are 
observed. 



• Modelling technique: the outputs of the modelling 
technique used for the analysis presented in this section, as 
described in detail in Paper I, are ranges of age, redden- 
ing and percentage contribution of both IYSPs and VYSPs. 
Therefore, part of the scatter in the YSP properties could 
be due to the uncertainties inherent in the modelling tech- 
nique; such uncertainties could hide underlying trends and 
patterns in behaviour. 

• Other variables: there are other variables such as, 
for example, the gas contents of the parent galaxies or the 
geometry of the merger, that could have a more important 
impact on the properties of the VYSPs in ULIRGs than the 
spectral classification, MFIR colours, MFIR luminosities, or 
the interaction class. 

Interestingly, this apparent lack of correlations between 
the various properties of ULIRGs has been found by other 
authors. Veilleux et al. (1999) found no correlation be- 
tween the emission line reddenings and other properties of 
ULIRGs, such as the MFIR colour ratio /25//60, for th e 108 
objects in their sample. Moreover, iFarrah et al.l l|200ll ) car- 
ried out an HST-WFPC2 V-band imaging study of a sam- 
ple of 23 ULIRGs and found no correlation between the IR- 
luminosity and the morphologies of the objects in their sam- 
ple. 



3 THE MASSES OF THE YSPS 

Table 4 in Paper I presents estimates for the total stellar 
masses associated with the different stellar components for 
the objects in the ES. Figure 2] shows the distribution of 
the total stellar masses presented in that table (account- 
ing for all VYSP, IYSP and OSP if present). Note that, 
in the cases of IRAS 14394+5332, IRAS 17028+5817 and 
IRAS 23327+2913, the total stellar mass was estimated for 
the two nuclei individually and, for the purpose of the fig- 
ure, the masses of the two nuclei were added. The figure 
shows that only 6 objects in our sample — IRAS 00188-0856, 
IRAS 13451+132, IRAS 14348-1447, IRAS 17179+5444, 
IRAS 23233+2819 and IRAS 23327+2913 — have stellar 
masses M st ciiar > 1 x 10 11 Mq. Although two of these ob- 
jects — IRAS 1345+132, IRAS 17179+5444 — are spectro- 
scopically classified as Sy2 galaxies in the optical, overall 
there is no clear-cut distinction between Seyfert and non- 
Seyfert ULIRGs in terms of their stellar mass; many Seyfert 
ULIRGs have stellar masses that are relatively modest - 
at or below the median for the sample as a whole. Thus 
our results provide no clear support for the idea that AGN 
activity in ULIRGs is more likely to be associated with rel- 
atively early-type galaxies with massive bulges, as appears 
to be the case in large samples of galaxies drawn from the 
Sloan Digital Sky Su rvey (SDSS) jkauffmann et all 12004 
iHeckman et all |2004 ). 

In order to compare the results found here with other 
studies, Table [6] shows the total mass values in units of m*, 
defined as the break in the measured mass function obtained 
from a large sample of 17,173 galaxies drawn from a com- 
bination of the Two Micron All Sky survey (2MASS) Ex- 
tended Source catalogue an d the 2dF galaxy redshift sur- 
vey (m. = 1.4 xlO 11 Mm: ICole et all l200ll , adapted to 
our cosmology). For the ES sample as a whole, the stellar 
masses are in the range 0.1 m» < M stc iiar < 3.6 m*, with a 
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Object Name 


Mgtellar 


Mi 

dyn 


Mi 

dyn 


M j 

dyn 


IRAS 


(m*) 


(m») 


(m«) 


(m*) 


00091— 0738 


0.33 








Oni 88-0856 


0.72 








01004-2237 


0.11 






0.07 


085724-3915 


0.20 




0.13 




1 01 QD-l-1 399 


0.28 






0.93 


1 04944-4424 


0.37 








12072-0444 


0.26 








1211 24-0305 


0.52 




0.70 


0.40 


1 25404-5708 










13305-1739 


0.28 








1 34284-5608 


0.53 


2.0 


0.6 - 1.1 


0.88 


134514-1232 


3.60 








1 35394-2920 


0.24 








140604-2919 


0.40 








14252-1550 


0.38 








14348-1447 


0.77 


1.0 


1.1 


1.50 


143944-5332 


0.44 








15130-1958 


0.28 






0.51 


1 52064-3342 


0.60 




0.3 




153274-2340 


0.25 


0.3 


0.3 


0.26 


161564-0146 


0.22 








1 R4744-3430 


0.61 








1 64874-5447 


0.67 








170284-5817 


0.40 








170444-6720 


0.19 








171794-5444 


1.24 








20414-1651 


0.25 






0.52 


21208-0519 


0.68 








22491-1808 


0.20 






0.68 


23060+0505 


0.49 








232334-2817 


1.64 








232344-0946 


0.34 








23327+2913 


3.37 








23389+0303 


0.15 









Table 6. Estimated total masses associated with the stellar populations (YSPs, including VYSPs and IYSPs, and OSPs if present) 
detected at optical wavelengths, expressed in units of m* (m» = 1.4 xlO 11 Mq). Due to the powerful AGN emission, it was not possible 
to estimate the total mass for the three objects in our sample classified as Syl galaxies (IRAS 12540-5708, IRAS 15462-0405 and IRAS 
21219-1757). 

1 Dynamical mass from lTacconi et alj 

2 Dynamical mass from Colina ct al. ( 20 05|). 

3 Dynamical mass from iDasvra et ail J2006al lbh . 



mean value of 0.6 m* and a median of 0.4 m». Aproximately 
80% of the objects in the ES (27) have total stellar masses 
M stc n ar < 1 xlO 11 M Q (0.71 m»). Also shown in the table 
are the dynamical mass estimates (inc l uding stars, gas and 
dar k matter) fromlricconi et all (|2002o , IColina et ail (|2005l ) 
and IDasvra et aL I i|2006al lbl) for the objects in common with 
the ES. The values presented in the table are for the entire 
systems, including the extended regions of the galaxies and 
adding the dynamical masses of the two nuclei for the dou- 
ble nucle us systems. To measu re th e dynamical mass of t he 
systems, iTacconi et all (|2002h and IDasvra et all (|2006al |bh 
used near-IR CO and Si I abs orption featur e s, trac ing the 
stars in the galaxies, although IDasvra et al.l (|2006al lbh also 
used the [Fe II] forbidden emission line, which i s a tra cer of 
the warm, ionized gas. In contrast. IColina et alj (|2005l ) used 
the optical Ha emission line, which also traces the warm, 
ionized gas. The main uncertainties associated with these 
techniques are: (1) the internal extinction; (2) the light con- 
tribution from AGN and/or compact, nuclear starburst; and 



(3) the assumption of v irialization (for a detailed discussion 
see IColina et all 120051 ). A comparison between our mod- 
elling results and those of the kinematical studies mentioned 
above shows that: 

- typically the stellar masses are within a factor of ~2 of 
the dynamical masses, which is remarkable given the uncer- 
tainties in both spectroscopic and dynamical techniques; 

- the dynamical masses tend to be larger than the stellar 
masses in most cases. However this is not surprising given 
that the dynamical masses account for all the mass (includ- 
ing stars, gas and dark matter), whereas the spectroscopic 
masses only account for the stars. 



The results presented here reinforce t he idea that 
ULIRGs are sub-m* or systems llGenzel et al 



l200ll:lTacconi et all. 120021 : IColina et all 120051 : IDasvra et al 



l2006al lbl). They also provide evidence that the optically vis- 
ible stellar populations dominate the stellar masses of the 
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Figure 4. Histogram showing the distribution of the total stel- 
lar mass values presented in Table [6] for the objects in the ES. 
The figure shows that only 6 of the objects have stellar masses 

Millar > 1 X 10 11 M . 



systems; the mass contributions of any stellar populations 
entirely hidden by dust must be relatively minor. 



4 BOLOMETRIC LUMINOSITIES: HIDDEN 
STARBURST COMPONENTS 

As well as the total stellar masses, the estimated values for 
the bolometric luminosities associated with the stellar pop- 
ulations detected at optical wavelengths are also presented 
in Table 4, of Paper I. For the ES, we find a mean value of 
L bo i = 0.66 x 10 12 L and a median of L bo i = 0.53 x 10 12 
Lq . Assuming that most of the optical light is absorbed and 
re-processed by dust, it is possible to compare these values 
with the mid- to far-IR luminosities of the soures. Figure 
[5] presents a histogram showing the estimated bolometric 
luminosities of the YSPs detected in the optical expressed 
as percentages of the mid- to far-IR luminosities of the 33 
sources in the ES considered here (all but the three Syl, for 
which it was not possible to estimate Lboi). We find that for 
48% of the objects (16 of 33) the bolometric luminosities of 
the YSPs detected at optical wavelengths represent a large 
fraction ( > 50%) of the mid- to far-IR luminosities of the 
sources. 

However, note that 6 of the 16 objects in which 
Lboi(YSP) represents a large fraction of the Lir are spec- 
troscopically classified as Sy2 galaxies in the optical. Due to 
potential AGN contamination, the contribution of the stellar 
populations to the optical light, and therefore the bolomet- 
ric luminosity associated with them, is less well constrained 
for these objects than for the other ULIRGs. Hence, it is 
possible that the true values of the stellar bolometric lumi- 
nosity presented in Paper I are significantly smaller. Exclud- 
ing these 6 objects from the group of 16, we find that in 23 
of the remaining 33 objects in the ES (70%) the bolometric 
luminosities of the YSPs represent only a modest fraction of 
mid- to far-IR luminosities (L yap /Li r <0.5). On the other 
hand, there remain 10 objects in the ES (~30%), with no 



50 100 

YSP proportion of infrared luminosity (%) 

Figure 5. Histogram showing the estimated bolometric luminosi- 
ties of the YSPs detected in the optical expressed as percentages 
of the mid- to far-IR luminosities of the 33 sources in the ES 
considered here. 

evidence of AGN activity, for which we find that the YSP 
luminosities represent a large proportion of their infrared 
luminosities. 

At this stage, it is important to mention that there are 
two major uncertainties associated to the results presented 
in this section: first, uncertainties in the selection of a par- 
ticular model that is then applied to the whole galaxy; sec- 
ond, the assumption that all the stellar light is actually ab- 
sorbed by, and heats the dust that radiates at mid- to far-IR 
wavelengths. In the former case, our model selection and as- 
sumptions are unlikely be substantially in error, since the 
stellar masses derived from the same models compare well 
with the dynamical masses. Of more concern is the second 
assumption: that all the light from the stellar populations 
heats the dust. In at least some cases, the reddening of the 
VYSP in the nuclear aperture is fairly modest (see paper 
I for details) and the IYSP reddening is constrained to be 
E(B-V) of 0.4 or less. In the extended regions not covered by 
the slit (but included in our bolometric luminosity estimates 
because of the aperture correction of the flux), the redden- 
ing is likely to be even less. Therefore it is possible that a 
substantial fraction of the starlight from the stellar popula- 
tions detected in the optical escapes each galaxy, and that 
the stellar bolometric luminosities presented in Table 10 in 
paper I and Figure [5] in this paper represent over-estimates 
of the ability of the starlight to heat the dust. 



5 THE NATURE OF THE PROGENITOR 
GALAXIES 

We discuss in this section the extent to which the stellar 
populations detected in the optical have been formed dur- 
ing the merger, and also the nature of the progenitor galax- 
ies involved in the triggering events. The aim is to answer 
two simple, but important, questions. First, are the YSPs 
detected in optical observations of ULIRGS all formed dur- 
ing the merger event? Second, can we use studies similar to 
that presented here in order to determine the star formation 

histories in maj or galaxy mergers? 

As outlinedfRodrigue z Zaurfn et al.ll|2008h . the contam- 
ination by YSPs present in the disks of the galaxies prior to 
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the merger is potentially an important issue. We investigate 
this by comparing the VYSPs masses of late-type spirals, 
likely progenitors of ULIRGs, with those of the IYSPs in 
the ULIRGs in our ES. The reason to focus on the VYSPs 
in spirals is that, if the star formation is truncated dur- 
ing the merger, these are the stellar populations that may 
evolve into the IYSPs observed today. For the ULIRGs in 
the ES, we obtain IYSP masses in the range Miysp = 0.5 
50 x 10 10 M B , with a mean value of 6.0 x 10 10 M and 

comparison, 



a median of 3.4 x 10 10 M ( see paper 10 In 



iRodrfguez Zaurfn et all (120081 ) estimated that the VYSP of 
typical late-type (Sc) spiral galaxies have VYSP masses in 
the range Mvysp(Sc) = 1.2 - 8.5 x 10 9 M Q . Therefore it is 
clear that, assuming an equal mass merger, both of the pro- 
genitor spiral galaxies must be at the upper end of the mass 
range in order to explain the IYSPs in ULIRGs in terms of 
VYSP from the captured disks of the progenitor galaxies; 
clearly the ULIRGs with the more massive IYSPs would be 
difficult to explain in this way. 

Alternatively, we can consider the nature of the pro- 
genitor galaxies that would be required to produce the 
YSP (IYSP and VYSP) and gas masses of ULIRGs 
by gas-rich mergers, with the YSP formed entirely by 
merger-induced star formation. The total gas masses M(H2 
+ Hi) of ULIRGs are typi c ally of order a few times 
10 10 M m jMirabel fc Sanders!. 1 19881 ; ISolomon et ail . 1 19971 ; 
lEvans et all 120021 ). Adding these gas masses to the total 
masses of the IYSP and VYSP (from Table 4 in paper I) 
that we are assuming have formed in the mergers, we find 
a total mass M{gas + YSP) > 5 x 10 10 M in many of 
the objects in our sample. Given that the merger-induced 
star formation is unlikely to be more than 50% efficient, the 
progenitor galaxies are together required to have an even 
larger gas mass. In comparison, typical late type (Sc) spiral 
galaxies have total gas masses of M(H2 + Hi) ~ 1.4 x 10 
Mq (see IRodrfguez Zaurfn et all [2008. for details). There- 
fore the progenitor galaxies are required to be at the upper 
end of the mass range for late-type spirals to explain the 
total YSP and gas masses of the ULIRGs. Note that, even if 
we consider only the gas masses, we require two "average" 
late-type spiral galaxies to merge in order to produce the 
typical gas masses measured in ULIRGs, and this is without 
accounting for any losses of gas due to outflows and star 
formation in the merger. 

Finally we note that 33% of the objects in the ES have 
VYSPs with masses Mvysp > 10 10 M (see paper I). Since 
these VYSP must be formed during the merger, and the star 
formation efficiency is likely to be <50%, the total amount 
of gas required in the progenitor galaxies is much larger; 
again a merger of two late-type spirals of average or larger 
mass would be required to produce the VYSP alone in these 
systems. 

Overall, regardless of whether the IYSP are formed dur- 
ing the mergers or merely captured from the progenitor 
disks, it is clear that the progenitor galaxies are required 
to be massive g as-rich spiral galaxies. Usi ng the compilation 
of HI masses in lRoberts fe Havnesl (| 19941 ). we find that, for 



the cases in which the estimated total progenitor gas mass 
is M(gas + YSP) > 5 x 10 10 M , both progenitor galaxies 
are required to be in upper 25% of the gas mass range for 
spiralqj. 

At this stage, it is important to add a caveat about 
the nature of the progenitor galaxies. It is commonly ac- 
cepted that the ULIRG phenomenon is associated with gas- 
rich mergers in which both galaxies are spiral galaxies. This 
is supported by the mass arguments presented above, as well 
as the lack of evidence for a major contribution from an OSP 
in most of the ULIRG sample. However, in the cases of one 
of the nuclear apertures in each of the double nucleus sys- 
tems IRAS 21208-0519 and IRAS 23327+2913, adequate fits 
are only obtained with a combination of stellar populations 
that includes a dominant contribution of a 12.5 Gyr OSP, 
plus a VYSP. Overall, these results suggest that mergers 
in which at least one of the progenitors had an early-type 
morphology can also trigger the ULIRG phenomenon. 



6 COMPARISON WITH MERGER MODELS 

The range of different morphologies of the objects in the 
ES, from widely separated double nuclei systems (e.g. IRAS 
14394+5332) to single nucleus galaxies with no strong signs 
of tidal structures (e.g. IRAS 15206-3342), suggests that 
they are in different stages of merger events. In this sec- 
tion, we give an overview of how our modelling results fit in 
with the merger simulation predictions. 

In general terms, the simulations predic t two 
epochs of starburst activity dMihos fc Hernguistl 1 19961 ; 
iBarnes fc Hernguistl Il996l ; ISpringel et al l . 120051 ) m major 
gas-rich mergers: the first occurring around or just after 
the first encounter, and a second more intense episode to- 
wards the end of the merger, ~0.5 - 1.5 Gyr after the first 
encounter, when the nuclei are close to coalescing (within 
5xl0 7 Myr of coalescence). However, both the time lag 
and the relative intensity of the peaks of starburst activ- 
ity during the merger event depend on several factors: the 
presence of bulges, feedback effects, gas content and or- 
bital geometry. For example, the presence of a bulge acts 
as a stabilizer of the gas against inflows and the forma- 
tion of bar structures, allowing strong er starburst activity 
towa r ds the end of the merger event (|Mihos fc Hernguistl 
ll996l ; lBarnes fc HernouTstlll996l ). On the other hand, AGN 
feedback effects (e.g. quasar-driven winds) disrupt the gas 
surrounding the black hole, acting against t he star formation 
activity at all stages (|Springel et al. . l2005l ). 

In Paper I we showed that it is possible to model all 
of the extracted spectra with a combination of two stellar 
populations, plus a power-law in some cases. Concentrat- 
ing again on comb III, the modelling results are, in general 
terms, consistent with the merger simulations. For example, 
it is possible that the IYSPs detected in most apertures are 
related to the first enhancement of the starburst activity, 
around or just after the time of the first encounter, and we 
are now seeing the objects at a later stage. 



4 It is important to emphasize that the upper limit of the range 
is found for the radio galaxy IRAS 13451+1232, which clearly 
represents an exceptional case among the ULIRGs in the ES 



5 Assuming an H2/H I mass ratio of 0.7 for spirals 
iYoung fc Knezekl . ll989ll . 
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We emphasize that, although the VYSPs are more sig- 
nificant in the nuclear regions of the galaxies, this compo- 
nent is detected at all locations sampled by the extended 
apertures for the overwhelming majority of the objects. We 
also find a trend for the VYSPs in the nuclear regions to 
be younger and more reddened than those of the extended 
regions. All of these results are consistent with the merger 
simulations, which predict that the star formation activity 
in the final stages of merger is concentrated towards the nu- 
clear regions of the merging galaxies, coinciding with higher 
concentrations of gas and dust, but it also occurs within the 
tidal structures formed during the interaction. Furthermore, 
as described in Section 2, when using Comb I, we find a clear 
cut-off in YSP ages at 100 Myr in the nuclear regions. This 
is in good agreement with the timescale (< 100 Myr) of 
the induced starburst activity predicted by the theoretical 
models for the final stages of a merger, as the nu clei coalesce 
<|Mihos fc Hernauistl 1 19961 ; ISpringel et ail . I2005T ). 

We note that ongoing star formation activity is also 
detected for all objects with relatively wide nuclear sep- 
arations: IRAS 14394+53323 (54 kpc), IRAS 17028+5817 
(25 kpc), IRAS 21208-0519 (15 kpc) and IRAS 23327+2913 
(24 kpc). Three of these objects, IRAS 17028+5817, IRAS 
21208-0519 and IRAS 23327+2913, show no signs of AGN 
activity. Star formation rates of SFR > 100 Mq yr" 1 are re- 
quired in order to produce the i nfrared luminosit ies (Lir > 
10 12 Lq) in these latter objects IjKennicuttl 1 1998 ) . As noted 
above, such high star formation rates are only predicted 
by the merger simulations in two stages during the merger 
event: 

• at, or after, the first encounter, if the parent galaxies 
have insignificant bulges and feedback effects are not impor- 
tant, or 

• at the end of the merger event, close (within ~50 Myr) 
to the coalescence of the nuclei. 

Estimates of the dynamical timescales for completion of the 
mergers ranges from >50 Myr in the case of IRAS 21208- 
0519, to >80 Myr in the cases of IRAS 17028+5817 and 
IRAS 23327+29133 Therefore, if the star formation activ- 
ity is associated with the final starburst that occurs around 
the time of the coalescence of the nuclei, then this starburst 
must be relatively long-lived in these three wide separation 
systems. Alternatively, the enhanced star formation activity 
in these systems may be associated with the (relatively ex- 
tended) period of star formation that follows as the galaxies 
move apart the first encounter of the two nuclei. However, 
for the latter case, the existing merger simulations do not 
predict the large star formation rates (> 100 Mq yr" 1 ) for 
the extended period after the first encounter, that would 
be required for these systems to be observed as ULIRGs. 
Most probably, these wide nuclear separation ULIRGs are 



6 Note that in the case of IRAS 14394+5332, the widest separa- 
tion system in the sample, the ULIRG phenomenon is not related 
with the interaction of the two, widely separated systems, but 
with the eastern source, which is itself a merger in its final stages. 

7 In making these dynamical timescale estimates we have as- 
sumed that the projected distances between the nuclei represent 
their true separation, and that the nuclei are moving radially to- 
wards each other with a relative radial velocity of 300 km s — 1 . 



systems in which the final starburst has been triggered rela- 
tively early in the period prior to coalescence, as the nuclei 
move together. 



7 EVOLUTIONARY SCENARIOS 

The presence of evolutionary links between ULIRGs, 
QSOs/radio galaxies and elliptical galaxies is commonly ac- 
cepted. However, the true nature of the links remain uncer- 
tain. While some authors propose the evolutionary sequence 
cool ULIRG — > warm ULIRG — > opt i cally selected QSO or 



Canalizo fc Stockton! 
20011 : iTacconi et al.l . 



radio galaxy (|Sanders et al. (1988a)l 
I2001T) . other authors dGenzel et al.l . 
2002) dismiss such an evolutionary sequence based on the 
location of ULIRGs in the fundamental plane, similar to 
that of the intermediate mass (~ L„) elliptical galaxies. We 
will compare here our study of a sample of 36 ULIRGs with 
studies of the stellar populations in samples of QSOs and 

rela ted objects. 

ICanalizo fc Stockton! (|2000al lbl. l200lT ) carried out a sim- 
ilar study to that presented here for their sample of 9 "tran- 
sition QSOs" — objects that may represent a transitionary 
stage between ULIRGs and QSOs. All of their sample ob- 
jects show clear signs of interaction, and five of them are 
classifed as major mergers. Furthermore, all but one have in- 
frared luminosities Lir > 10 12 Lq (i.e. they are ULIRGs). In 
order to study a possible link between these objects and the 
general populations of ULIRGs, Canalizo & Stockton inves- 
tigated the stellar populations within the host galaxies. For 
that purpose, they modelled the optical spectra extracted 
from a series of extended apertures using a combination of 
a 10 Gyr OSP plus YSPs with varying ages, which is similar 
to the Comb I used here. 

Figure|6]presents histogr ams comparing the di s tributi on 
of th e YSP ages derived by ICanalizo fc Stockton! (|2000al lbl. 
l200lf) for the transition QSOs, with those obtained for the 
ULIRGs in the ES sample using Comb I (labelled as LW- 
age in Section 3). For the transition QSOs, we hav e used the 
starburst peak ages (ICanalizo fc Stocktonl l200ll . column 5 
of their Table 3), defined as th e predominant starbu r st age 
found in the host galaxy (see ICanalizo fc Stocktonl Il997l . 
l2000al lbl for details), while for the ULIRGs we have used 
the average values of the LW-age, i.e. the same as in Figure 
[T] Note that, since the nuclear regions could not be sam- 
pled due to contamination by the bright quasar nuclei, the 
Canalizo & Stockton studies concentrate on the extended 
regions of the objects. Therefore, we have compared their 
results with those for the extended regions of the ULIRGs 
in the ES. Figure [6] also shows the results obtained in this 
study for the extended apertures of the Sy-ULIRGs (the 
"oldest" ages in this figure are obtained for the atypical ob- 
ject PKS1345+12). 

As described in Section 2, we do not find any signifi- 
cant differences between the ages obtained for the HII- and 
LINER-ULIRGs and those of the Sy-ULIRGs in the ES. The 
same conclusions are reached if we compare the transition 
QSOs with the Sy-ULIRGs in the ES sample, or with the 
sample of ULIRGs as a whole. If there is an evolutionary 
link between ULIRGs as a class and the transition QSOs, 
one would expect the latter to have, in general, older YSP 
ages than ULIRGs, although some overlap is expected. How- 
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Figure 6. The young stellar age distribution for the Canalizo & Stockton (2001) transition QSOs (left) and the extended regions of the 
ULIRGs in the ES (right), with Seyfert ULIRGs highlighted. These results are based on modelling combination I. It is clear that the the 
transition QSOs, Seyfert ULIRGs, and non-Seyfert ULIRGs overlap in terms of their distributions of YSP ages in these diagrams. 



ever, it is clear from the figure that the ages obtained for the 
extended apertures in the transition QSOs are very similar 
to those of the extended apertures of the ULIRGs in the ES. 
A first interpretation of this result is that there is no clear 
evolutionary link between non-Seyfert ULIRGs and transi- 
tion QSOs. However, as mentioned in section 2, it is possible 
that the timescale for the transition from ULIRGs to QSOs 
is smaller than the timescale associated with the ULIRG 
phenomenon. 

Therefore, our study does not provide any clear-cut ev- 
idence for an evolutionary link between ULIRGs and the 
QSOs. Certainly, it is unlikely that the AGN in the Seyfert 
ULIRGs in our sample have been triggered a substantial 
period (> 100 Myr) after the main merger-induced star- 
burst, as appear s to be the case i n som e radio galaxies 
jTadhunter et all . 120051 ; IWills et al l . I2008T ) 



8 COMPARISON WITH STUDIES OF HIGH-Z 
OBJECTS 

Star-forming galaxies have been detected in large numbers 
at high redshifts (z > 1). These galaxies are selected on the 
basis of t heir strong UV em ission (Lyman-break galaxies - 
LBGs: see lGiavaliscol . [200l. for review), high sub-mill i meter 
luminosities (sub-mm galaxies - SMGs, ISmail et all 20021 ; 
iBlain et all |2002| ; IChapman et~ai] , 120031 ; iPope et all l2005j ) 
mid-IR flux densities (2 4/xm Spitzer- s electe d galaxiejf , 



Papovich & et al , 
20071 ; ISaiina et all 



2004 



Caputi et all . 
20071 ; iFarrah et 311 12008) 



20061 ; lYan et al 

or red near- 

IR col ours: (J — - fQvega ^ 2.3 (distant red galaxies - 
DRGs. iFranx et all 120051 ). A large fraction show morpho- 
logical features sugg esting that major mergers are common 
among such objects (iGiavalisc o. 2002; C hapman et allfeOQol ; 
lErb et all 120031 . |2004| ; IConselice et all 120051 ). In addition. 



8 Note that within the so called 24^m Spitzer-selected galaxies 
class there will be galaxies classified as LBGs, SMGs and DRGs. 
However, not all LBGs and/or DRGs will be detected at 24/xm 
and therefore, we refer to the 24/^m Spitzer-selected objects as a 
class of galaxies with their own properties. 



Star forming 


Medlar 


references 


galaxy type 


10 10 


Mq 




Lyman-break 


0.2 - 


6.0 


Papovich et al. (2001) 


Sub-mm 


12 - 


25 
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Table 7. The stellar masses of the different types of high-z star 
forming galaxies found in different studies. The stellar masses for 
the ULIRGs found in this study are also shown in the table for 
comparisson. Note that the upper limit of 50xl0 10 M Q is fouind 
for the radio galaxy PKS1345+12, wich represents an exceptional 
case among the ULIRGs in our sample. All results obtained for 
the high-z star forming galaxies are based on the modelling of 
photometric points from the optical to the mid-IR. 



iFraver et all (2004), based on deep near-IR imaging, con- 
cluded that the colours, luminosities, morphologies and sizes 
of the SMGs in their sample are all consistent with the prop- 
erties of local ULIRGs. Furt hermore, on the bas is of their 
CO millimeter observations, iTacconi et al. I (|2006h conclude 
that SMGs appear to be the "scaled-up" version of local 
ULIRGs. Therefore, it is interesting to compare the results 
of this study with those obtained for high-z star forming 
galaxies. 

Table [7] shows a compilations of the results obtained 
from various studies of the stellar masses of high-z star form- 
ing galaxies. We find that the masses of the local ULIRGs 
in our sample are, in general, comp arable with or slightly 
higher than, those of the LBGs in the lPapovich et al] (1200 ll ) 
sample. In the cases of the SM Gs and the DRGs, the resu lts 
found bv lTacconi et al.. 20081 and lPapovich et al. (2006)1 are 
consistent with the total stellar mass estimated in this work 
for local ULIRGs, and reinforce the idea that the latter ob- 
jects are the low redshift analogues of the SMGs and DRGs 
found at higher redshifts. On the other hand, we find that 
the estimated stellar masses for the ULIRGs in our sam- 
ple are significantly smaller (with the exception of the ra- 
dio galaxy PKS12345+12) than many of the Spitzer galax- 
ies. This resuts is c onsistent w i th tho s e of the recent stud - 
ies of, for example, lYan et al l (|2007h . ISaiina et all (|2007h . 
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ISaiina et all (|2008T l. iRigbv et all (|2008h and iFarrah et all 
(2008) among others, suggesting that the Spitzer galaxies 
are either maximum starburst, i.e. scaled-up versions of lo- 
cal ULIRGs, or indeed, represent a separate class of objects 
from local ULIRGs. 

Finally, note that the results obtained for the stellar 
populations of high-z star forming galaxies are based on 
fits to relatively few photometric points, albeit with a wide 
spectral range. The modelling work detailed in Paper I em- 
phasized the difficulty in obtaining a unique solution when 
using combinations of different stellar populations to model 
spectra, especially on the basis of fits to the SED alone (i.e. 
not examining the detailed absorption features). We have 
also shown that some results, and therefore the conclusions 
based on them, can change significantly with different mod- 
elin g assumptions. Further more, the SED fits for DRGs in 
the IPapovich et al. (2006)1 sample are relatively poor (x 2 > 
2) in a large fraction ( > 50%) of the sources. Such poor fits 
would not be accepted as viable in our study of ULIRGs. 
Clearly, particular care must be taken when interpreting the 
results obtained on the basis of modelling a few photomet- 
ric points using combinations of stellar templates, which is 
often the technique used in high-z studies. Studies similar to 
the one presented here are vital for increase our undertand- 
ing, not only of the properties of the stellar populations in 
high-z star forming galaxies, but also the nature of the links 
between these objects and the local ULIRGs. 



9 SUMMARY AND CONCLUSIONS 

In this paper we have presented a detailed analysis of the 
modelling results presented for the CS and the ES samples in 
Paper I. In addition, we have also discussed these results in 
a more general context, comparing them with other studies 
of ULIRGs. The conclusions can be summarized as follows. 

• Age, reddening and percentage contribution: we 

find that the YSPs are more significant in the nuclear than in 
the extended regions. The statistical analysis presented here 
suggests that the YSPs located in the nuclear regions tend 
to be younger and more reddened, although further analysis 
is required to confirm this result. All of these results are 
consistent with the merger simulations. 

• Correlations: we find no evidence for correlations be- 
tween the properties of the stellar populations and other 
properties of ULIRGs. This can be explained in terms of 
a selection effect, i.e. since the objects in the CS and the 
ES samples are ULIRGs, with high IR-luminosities and star 
formation rates, it is likely that we are observing them at, 
or after, the first encounter, or at the end of the merger 
event, when the nuclei are very close to coalescing. There- 
fore, they all have similar properties. A second possibility 
to explain the apparent lack of correlations is the scatter in 
the properties of the YSPs due to uncertainties in the mod- 
elling technique. A third possible explanation is that other 
variables, such as gas content of the parent galaxies and the 
geometry of the collision, may have a more important im- 
pact on the stellar properties of ULIRGs than the optical 
spectral classification, interaction class or luminosity. 

• Merger simulations: in general our modelling results 
are consistent with the merger simulations, with the IYSPs 



in most systems plausibly associated with the first enhance- 
ment of activity following the first encounter of the nuclei of 
the merging galaxies, and the VYSPs likely related to the 
final enhancement as the nuclei merge together. 

• 100 Myr cut-off: when using Comb I, we find a clear 
cut-off in luminosity weighted YSP ages at 100 Myr for 
the nuclear apertures, consistent with the timescale of the 
merger-induced starburst activity predicted by the merger 
simulations in the final stages of the merger event, as the 
nuclei coalesce. This result supports the idea that, despite 
the varying morphologies, we are observing most of the ob- 
jects close to the peak of the star formation activity in the 
final stages of the merger event. 

• The masses of the parent galaxies: both progenitor 
galaxies in the trigger events must be in the upper 25% of 
the mass range of late-type spirals for most of the ULIRGs 
in the ES. 

• The morphologies of the parent galaxies: although 
in most cases the results are consistent with the idea that 
the progenitors are late-type spirals (albeit unusually mas- 
sive late-type spirals), in 3 cases (IRAS 13451+1232, IRAS 
21208-0519 and IRAS 23327+2913) there is clear evidence 
that at least one of the progenitors is an early-type galaxy. 
This demonstrates that ULIRGs can potentially be triggered 
by mergers between galaxies with a range of galaxy types. 

• Stellar mass content: our modelling results suggest 
that most ULIRGs are sub-m*, or m* systems (m* = 1.4 
xlO 11 Mq). Such masses support the idea that these sys- 
tems will eventually evolve into intermediate mass elliptical 
galaxies. The masses obtained from the modelling results 
are generally consistent with the dynamical mass estimates 
(when available) within a factor of ~ 2. This result shows 
that the stellar populations detected in the optical dominate 
the stellar masses of the galaxies. 

• Bolometric luminosities: our results suggest that the 
stellar popultions detected at optical wavelengths make a 
significant contribution to the mid- to far-IR luminosities 
of many of sources. However, due to uncertainties related 
to the combination of stellar populations assumed, and to 
the assumption that all the stellar emission is absorbed by, 
and heats, the dust, it is possible that the values presented 
in Table 10 in paper I and Figure [5] in this paper are over- 
estimates. We conclude that while the visible stellar popu- 
lations can make a significant contribution to the heating of 
the dust, it is unlikely that they dominate the heating of the 
dust in most cases. 

• The evolution of ULIRGs: the stellar masses of the 
ULIRGs in the ES sample are consistent with the idea that 
ULIRGs evolve into intermediate mass elliptical galaxies. On 
the other hand, the data do not provide any clear evidence 
to support/refute the evolutionary sequence: cool ULIRG 
— > warm ULIRG — + optically selected QSO, since the stellar 
populations of the cold/warm, Sy/non-Sy, QSO/non-QSO 
ULIRGs are similar (perhaps the timescale of this transition 
is shorter than the timescale of the major merger-induced 
starburst). However, it remains possible that some ULIRGs 
evolve into radio galaxies (although not all radio galaxies 
can have this origin). 

• High-z counterparts: many of the properties of the 
ULIRGs in our sample, such as morphology, stellar mass 
content, star formation histories etc., are consistent with the 
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idea that ULIRGs are the local counterparts of the sub-mm 
galaxies (SMGs) and distant-red galaxies (DRGs) . 
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